Detailed study of a transition point 
in the Veneziano-Wosiek model of Planar Quantum Mechanics 

Piotr Korcyl 

M. Smoluchowski Institute of Physics, J agiellonian University 
Reymonta 4, 30-059 Krakow, Poland 

Abstract 

Following a model recently investigated by Veneziano and Wosiek we briefly introduce Planar Quantum 
Mechanics (PQM). Then, we present high precision numerical results in the sectors with two and three 
fermions. We confirm, that the transition point in the 't Hooft's coupling constant A in these sectors occurs 
at Ac = 1, as was expected. 
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The exact numerical spectra of a branch of re- 
duced supersymmetric theories can be calculated in 
a cut Fock basis by a method proposed recently by 
Wosiek [T] . In a series of papers [H [3] he studied some 
models of Supersymmetric Yang-Mills Quantum Me- 
chanics (SYMQM). These systems result from a dimen- 
sional reduction of the full dimensional (D = d + 1) 
supersymmetric Yang-Mills quantum field theories to 
a single point in space (0 + 1). The method pro- 
vided a good understanding of the D — 2 and D = 
4, N = 2 spectra [U [3]. The goal of such analysis is 
to reach, on one hand the system with D = 4, = 3 
which could shade some light on the real QCD, and 
on the other, the D = 10, SU(iV — > oo) model, which 
is conjectured to be in relation with the M-theory [4j. 
The latter case needs to incorporate into the scheme 
the large N limit [5], which is expected to provide a 
remarkable simplification. It should reduce consider- 
ably the number of basis vectors to be taken into ac- 
count and allow to disregard all non-planar contribu- 
tions. The above was investigated by Veneziano and 
Wosiek in [51 [71 [5] on a simple supersymmetric model 
belonging to the class of Planar Quantum Mechanics 
(PQM). 

In this paper we present high precision numerical re- 
sults on the model from Ref[6] in the sectors with two 
and three fermions. We investigate the transition point 
in the 't Hooft's coupling A, and show that it indeed 
takes place at Ac = 1. Therefore we give a numeri- 
cal confirmation of the analytical results obtained by 
Beccaria in [5] . 

The paper is constructed as follows. We start with 
an introductory part to the PQM, then we present the 
studied model, and finally we discuss our results. 

Supersymmetric Yang-Mills 
Quantum Mechanics 

We will introduce now the systems called SYMQM 
in the Hamiltonian formulation. Let's consider a 
quantum mechanical system with N'^ bosonic and iV^ 
fermionic degrees of freedom. As was already men- 
tioned, it can be regarded as the remainder after a 
dimensional reduction of supersymmetric field theory 
with U(iV) gauge symmetry to one point in space. Dur- 
ing such a procedure, a local gauge symmetry becomes 



a global one. Thus, our system should be invariant un- 
der a global V{N) rotation. Let T°- be the generators 
of the V{N) group in the fundamental representation, 
thus, they are N x N matrices. We introduce bosonic 
and fermionic matrix-valued annihilation and creation 
operators 

a.,^a%, al^a^%, (1) 

h = f''T^^, 4 = /^'^' (2) 

where the sum over b — 1,...,N'^ is assumed, and 
J, J = 1, . . . , iV. The invariance of the system is assured 
by taking the Hamiltonian as a trace of a polynomial 
of the above operators. The creation and annihilation 
operators satisfy the following commutation and anti- 
commutation relations 

^ Stl^jk, (3) 

{hJli}^SuS,k. (4) 

The Fock basis is composed of eigenstates of the 
occupation number operators, B = Tr{a^a) and 
F = Tr{pf), which are explicitly U(iV) -invariant. 
The construction of the basis starts from the Fock 
vacuum denoted by |0 >. We act on the latter with in- 
variant 'bricks' i.e. creation operators contracted with 
U(A^) invariant tensors. For the U(2) group we have 
two such tensors, the Sij and eijk- The basis states are 
obtained by an action of any combination of powers of 
theses bricks. 

The cut-off" method 

As it is impossible to handle infinite matrices on a 
PC, one needs to cut them somehow. The most in- 
tuitive way to do this is to introduce some integer, 
Bmax, and to keep only those basis states for which 
the total bosonic occupation number does not exceed 
Bmax- The method of obtaining the spectrum simply 
consists of calculation of the hamiltonian matrix ele- 
ments in such a cut Fock basis, and then of its numer- 
ical diagonalization. This should be done for several 
different cut-offs and a limit of infinite cut-off should 
be extrapolated in order to obtain a Bmax indepen- 
dent, thus physical, results. The difficulty of such a 
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program is hidden in the number of basis states grow- 
ing exponentially with increasing Bmax and A'^. Up to 
now, calculations have been made up to U(4) [TU] , 

Large N limit and Planar Quantum Mechanics 

The difficulties described in the preceding paragraph 
largely disappear in the 't Hooft limit, TV oo, 
g^N = const, where g is a coupling constant present 
in the system. The zeroth order approximation in the 
expansion consists in retaining only those contribu- 
tions which correspond to planar graphs. It appears 
that this can be done already on the level of the Fock 
basis [6] . The main contribution will be given by basis 
states obtained by action of a single trace brick. There- 
fore, in the purely bosonic sector, F — 0, one needs to 
consider only basis states of the form 



|0,n >= 



1 



^fo,n 



Tr[(at)"]|0>, 



(5) 



which are labeled by one integer, n. We can calculate 
an explicit expression for the normalization constant 
A/o.Ti [6j. The sector F — 1 contains one fermion and 
the basis states are given by 



11,71 > = 



1 



Tr[{a^rp]\0> 



(6) 



With increasing fcrmionic occupation number F, 
things get complicated, because one has to use several 
integers to label basis states. For example, if = 2 
we need two integers, ni and n2- The basis state is 
therefore obtained by the action of a trace [7] 



|2,ni,n2 > = 



2,711 -^2 



-Tr[{a^r^fHa^rf]\0> . (7) 



Due to the cyclicity of the trace, we only need to 
deal with states with ni < n2- Moreover, if ni ~ 
7X2 = m, the anticommutation of fermionic creation 
operators and the cyclicity of the trace imply that 
|2, m,m >= 0. So, the basis is composed of states for 
which ni < n2. Thus, for a given cut-off B^ax, we will 
have ^B^axiBmax - 1) states. 

In the case of three fermions, we need three integers 



to label a basis state 

\3,ni,n2,n3 >= (8) 
-rr[/t(at)"V^(a^)"V^(a^)"']|0 > . 



1 



•^2^111 ,n2 .Ti 

Again, we can arrange them so that ni < n2,n^. 

The large N limit of SYMQM systems in the ze- 
roth order approximation is called Planar Quantum 
Mechanics. 

Veneziano-Wosiek model 

The model considered in Ref[6l |7l [8j is given by the 
supersymmetry generators 

Q^Tr[fa\l+ga% = Tr[f\l + ga)a], (9) 

where g is the coupling constant. We can define the 
't Hooft's coupling constant as A = g'^N, where TV pa- 
rameterizes the gauge group \]{N). The Hamiltonian 
reads 

Hb = Tr[a)a + g{a^^a + a)a^)+g^a^^a\ (10) 
Hf - Tr[/t/ + g(/t/(at+a) + /t(at+a)/)(ll) 
+ g^f^afa^ + paa^f + pfa^a + pa^fa)]. 

It conserves the fermionic occupation number 
F = Tr[f^f], so we can analyze our model separately 
for each fixed F. The cases F = and F — 1 were 
described in Ref[5], whereas the sectors F = 2 and 
F = 3 in Ref[7j. We will concentrate here exclusively 
on these higher-fermion-number sectors. 

Following the rules of planar calculus, described in 
detail in Ref[5], one can calculate the matrix elements 
of the Hamiltonian in the sectors with two and three 
fermions. We just recall here the explicit results [7]. 

Two fermion sector 

We use the notation for the matrix element: 



,n2 ;mi ,m2 

= < 2,ni,n2\H\2,mi,m2 > 



Then 



H„ 



{ni +n2 + 2)(1 + A) - A(2 - (5„,,o + 2Sn^,n,+i) 



3 



-^^ni + l,n2;ni,n2 — -^ni ,712 ;ni + I,n2 ^ \/~X{tLi -\- 2) (13) 
-^Tii ,n2 + l;ni ,712 -^ni ,n2 ;ni ,n2 + l ^ V'^(7l2 4" 2) 



2xa - 



rii,n2:ni+l,n2 — 1 



1) 



(14) 



Three fermion sector 

Similarly, we denote the matrix element by 

,n2 ,713 'jl^l -,'^2 

= < 3, rti, n2, nsliJjS, mi, 1712,1123 > 

We have 

-^ni ,n2 -ns :ni ,n2 ,n3 — (-^^) 
(rii + 712 + «3)(1 + A) - A(3 - (5„,,o - (^na.O - (^ria.o), 



ni + l,n2,"3;"l,n2,n3 



Til, 712,713 ;ni + 1,712,713 



\/A(?ii +2)A, 

-^7ii,n2 + l,n3;ni,7i2,'n'3 — -^tii ,712 ,^^3;^! ,7i2 + l,' 



\/A(?i2 + 2)A, 



(16) 



H, 



ni,n2,n3 + l;ni,n2,n3 — -^"l ,"2 ,«3 ,"2 ,"3 + 1 

%/A(n3 + 2)A, 



AA, 

AA, 

AA, 
(17) 



-ffni+l,n2 — l,"3;ni,"2,n3 — -f^rii ,n2 ,"3 ;ni + l,«2 — l,n3 
-^ni ,n2 + l,n3 — l;ni ,n2 ,n3 — -^^ni ,n2 ,n3 ;ni ,n24■l,?^3 1 
-^ni — 1,712 ,n3 + l;ni ,712 ,n3 -^ni ,712 ,713 ;ni — 1,712 ,7i3 + l 

where A is defined by 

if for the right state J^i = 712 = ^3; 
\/3 if for the left state 111—122—113, 
1 otherwise. 

Previous investigations showed the existence of a 
transition point in the 't Hooft's coupling constant A 
at Ac = 1. On one hand, it appears as a critical slow 
down of the convergence of eigenenergies as a function 
of Bmax, and on the other, the spectrum becomes con- 
tinuous, whereas it was discrete away from Ac- It was 
possible to derive the existence of this transition point 
analytically in the sectors with none or one fermion. 
Numerical results strongly suggest that the transition 
also occurs at Ac = 1 in the higher-fermion-number 
sectors. The aim of the present paper is to confirm 
this by new high precision results from larger cut-offs 
calculations. 



High cut-ofF results 
and the transition point 

Our main goal here is to study in detail the location of 
the transition point Ac in the sectors with two or three 
fermions. Since the bases in these sectors are much big- 
ger than the ones in lower-fermion-number, one needs 
another tool for more quantitative analysis. We used 
ARPACK, a Fortran?? library for spare matrices, to 
diagonalize our Hamiltonian matrix. In this way we 
were able to reach cut-offs Bmax = 500 (110), respec- 
tively for F = 2 (3), corresponding to the sizes of basis 
up to 100000 vectors, compared to B,,nax = 40 (30) 
attained in Ref|?j. 

Two fermion sector 

We will find the transition point by examining the de- 
pendence of the energy of the ground state on the cou- 
pling constant A. Figure [T] shows this energy, which in 
the following we will call £'b,„„^(A), for a given cut-off 
Bmax and in some interval around A = 1. Suggestions 
from previous works are confirmed. Namely, for A > Ac 
the ground energy vanishes, and thus constitutes one 
of two SUSY ground states, which are present in this 
sector. For A < Ac, Eb^^^{X) is non-null and has a 
nontrivial dependence on A. The determination of the 
transition point is carried out by fitting to -E_b„„^ (A) a 
polynomial in A for several fixed Bmax- This polyno- 
mial is chosen to be positive for A < Aq and equal to 
zero at A = Aq 

w{X) = wi{X - Ao) + u;3(A - Aq)^ + W5{X - Xof. 

The fitted curves, together with the polynomial roots 
^o{Bmax), are shown in figure^ In order to obtain the 
value of the physical transition point Ac, we extrapolate 
^a{Bmax) to the limit Bmax 00. We do this by 
fitting two types of decreasing functions 



n 

7) 



Ac + b{B 

max J 

Ac + bexp{cBr, 



The resulting fits are shown in figure [5l whereas table [T] 
contains the values of fitted parameters. We can read 
off the infinite-cut-off limit Ac equal to 



Ac 



1.0061 ± 0.0005, 
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where the uncertainty is given by the difference be- 
tween the two values of Ac coming from the fits of the 
two functions. Finally, we also check the convergence 



Fitted function 


Obtained parameters 


Xo{Bmax) — 
Ac + b{BmaxY 


Ac = 1.0059 
b = 471.1 
c = -2.40 


Xo{B„iax) = 

Ac + b ex.p{cB„iax) 


Ac = 1.0064 
b = 0.060 
c = -0.022 



Table 1 : Numerical values of fitted parameters for the 
Xo{Bmax) dependence. 
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of the ground energy Eb„^^^ (A) at the suspected value 
of transition point A = Ac = 1.0. To this end, we cal- 
culate the extrapolation of the values Eb,^^^{\ = 1.0), 
obtained for some specific cut-offs, by fitting a function 

Eb„^^^ (A - 1.0) = Ec + biB^axf. 

The results for the fitted parameters are shown in 
table [21 and the curve is plotted in figure [6] We can 
conclude, that 

Ec = -1.08 10"^ ±4.5 10"^ 

where the error is given by the difference between the 
two results with highest cut-off. 

To summarize, in this paragraph we have showed, 
that in the sector with two fermions, the transition 
point occurs for Ac = 1.0 and that the energy of 
the ground state at the conjectured transition point 
A = 1.0 converges to zero. 



Fitted function 


Obtained 
parameters 


^^B„..(A = 1.0)^= 

Ec + h{BmaxY 


Ec = -1.08* 10"'^ 
6 62 
c= -2.73 



Table 2: Numerical values of fitted parameters for the 
^B^a^ci^ ~ 1-0) dependence. 



Figure 1: The dependence of the ground state of the 
sector with two fermions on the coupling constant A 
for different cut-offs B„iax- The highest dashed curve 
represents the results for the smallest Bmax = 100, 
whereas the lowest one corresponds to the highest 
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Figure 2: Polynomials fitted to the dependence of the 
energy of the ground state in sector with two fermions 
on the coupling A for different Bmax as well as their 
zero points. The highest dashed curve represents the 
results for the smallest B,nax — 100, whereas the lowest 
one corresponds to the highest B^ax = 400. 

Three fermion sector 

The analysis of the transition point in this sector fol- 
lows the lines of the preceding paragraph. Similarly, 
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Figure 3: The dependence of the ground state of the 
sector with three fermions on the coupling constant A 
for different cut-offs B^ax- The highest dashed curve 
represents the resuhs for the smallest B^ax = 45, 
whereas the lowest one corresponds to the highest 

Braax — 90- 

we will examine the dependence of the energy of the 
ground state, called £^b„„^(A), on the coupling A. Fig- 
ure [3] demonstrates the numerical curves for different 
cut-offs. The highest, dashed, curve represents calcula- 
tions for B„iax — 45, and the lowest one for Bmax = 90. 
We see that the convergence is very good for A away 
from Ac i.e. A < 0.90 and A > 1.05. The transition 
takes place for A between these values, and can be seen 
on this plot as a slow down of the numerical method. 
Let's denote, for each Bmax, the minimal energy of 
the ground state by Ejnin{Bmax) and its position by 
^min{Bmax)- The physical results, i.e. cut-off inde- 
pendent, are thus the limiting quantities and Ac 

such that Emin{Bmax) * Ef, and Xm,in{Bmax) * •^c 

Bmax oo. The results from sectors F = 0, 1 suggest 
that the energy of all states collapses to zero and we 
get a continuous spectrum at the speculated transition 
point Ac = 1. 

We determine \nin{Bmax) and Ermn{B,nax) by two 
methods. First of them consists in fitting a fourth order 
polynomial, 

w(A) = Wo -I- wiA + W2y? + WiX^, 




I ^ ^ ^ 1 

0.85 0.9 0.95 1 1.05 

A 



Figure 4: Polynomials fitted to the dependence of the 
energy of the ground state in sector with three fermions 
on the coupling A for different Bmax- The highest 
dashed curve corresponds to Bmax = 45 and the lowest 
one to Bmax = 90. The calculated minima are shown 
as well. 

to £'s,„„^(A) for each Bmax- Then, XminiBmax) and 

Emin {Bmax 

) are calculated analytically given the fitted 
parameters. The fitted curves are shown in figure [H 
together with the calculated minima. The second 
method uses the cubic spline to transform Eb^^^ (A) 
into a continuous curve. The approximated values of 
the minima are then found numerically by bracketing. 
The differences between the results coming from these 
two methods will be later used as an estimate of the 
uncertainty of the calculated quantities. To get Ec we 
extrapolate Emin{Bmax) to Bmax oo, and to this 
end, we fit a polynomial function 

Emin(.Bmax^ Eq ^" ^{Bmax^ - 

The obtained fit is shown in figure[7l whereas the values 
of parameters are presented in table [31 We thus have 

Ec = -0.00094 ±0.00021. 

In order to extrapolate XminiBmax) we fit three slowly 
growing functions: 

* Xmin{Bmax) ^ {Bmax) : 

* Xmin{Bmax) ^ {Bmax) ^ 
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Fitted function Obtained 

parameters 
ErmniB^ax) = ^ -0.00094 ± 0.00021 

Ec + b{Bma^,y 6 = 33.2 ±2.4 

c= -1.694± 0.019 



Table 3: Numerical values of fitted parameters for the 
ErmniBmax) dependence. 



Fitted functions 




Obtained 
parameters 






= 1.0155 ± 0.0053 


Ac + h{BjnaxY 


b = 


-0.474 ±0.045 




c = 


-0.480 ±0.043 


^mini^-B max') — 




1.01304 ±0.00038 


Ac + b{Bmax)~^^'^ 


b = 


-0.4964 ± 0.0033 


^mini^^max) 


Ac 


= 1.073 ±0.010 


Ac + blnic * Brnaxy^ 


6 = 


-0.478 ±0.084 




c 


= 0.80 ±0.28 



Table 4: Numerical values of fitted parameters for the 
>^min{.Bmax) dependencies. 



Discussion and conclusions 

In this paper we used high precision numerical re- 
sults in order to check the transition point in the 't 
Hooft's coupling constant A in the Veneziano-Wosiek 
model. We investigated the sectors with two and three 
fermions. By fitting some specific functions we extrap- 
olated from the numerical data the physical, i.e. cut-off 
independent, values of the transition point and ground 
energy at A = 1.0. We confirmed that in both sectors 
this transition point occurs nearly at A = Ac = 1.0, and 
that the ground energy at this value of coupling con- 
stant converges almost to zero. The uncertainty given 
with these results is not a true statistic error since it 
was not calculated from any statistical ensemble. It 
should be only interpreted as an indication of the real 
uncertainty. 
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Table |4] contains the obtained values of the fitted pa- 
rameters, and the fitted curves are shown in figure [S] 
Eventually, we can assume that the value of the con- 
stant coefficient Ac is equal to the mean of the values 
obtained from the three fits, and its error is the stan- 
dard deviation. Therefore 

Ac = 1.034 ±0.016. 

One also notes, that the general fit of a power function 
gave approximately the same results as the fit of the 
inverse of the square root. 

As a conclusion of this section we recapitulate our 
results for the sector with three fermions. Namely, we 
showed that the transition point occurs at Ac = 1.0 
and that at this value of coupling constant the ground 
energy converges to zero. 
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Figure 5: The fitted dependencies of \o{Bmax) on Bmax- 
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Figure 6: The fitted dependence of Eb^„^{X = 1-0) on Bmax- 
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Figure 7: The fitted dependence of Emin{Bmax) on Bmax- 
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Figure 8: The fitted dependencies of \mm{Bmax) on B^ax- 
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